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ABSTRACT 

We present an analytic method for rapidly forecasting the accuracy of gravitational potential recon- 
struction possible from measurement of radial peculiar velocities of every galaxy cluster with M > Mth 



in solid angle 9 and over redshift range 



< Z < Zn 



These radial velocities can be determined 



from measurement of the kinetic and thermal Sunyaev-Zeldovich effects. For a shallow survey with 
0.2 < z < 0.4, coincident with the SDSS photometric survey, one mode of the gravitational potential 
(on length scales > 60Mpc) can be reconstructed for every 8 cluster velocity determinations. Deeper 
surveys require measurement of more clusters per S/N > 1 mode. Accuracy is limited by the "under- 
sampling noise" due to our non-observation of the large fraction of mass that is not in galaxy clusters. 
Determining the gravitational potential will allow for detailed study of the relationship between galaxies 
and their surrounding large-scale density fields over a wide range of redshifts, and test the gravitational 
instability paradigm on very large scales. Observation of weak lensing by large-scale structure provides 
complementary information since lensing is sensitive to the tangential modes that do not affect the 
velocity. 

Subject headings: cosmology: theory - cosmology: observation - cosmology: weak lensing - cosmology: 
peculiar velocities - galaxies: formation - galaxies: evolution 



1. INTRODUCTION 

Although there are a variety of techniques for measuring 
the statistical properties of cosmological density fluctua- 
tions, we know of only three types of observations from 
which the density field itself may be reconstructed: weak 
lensing (Mellier, 1999; Bartelmann & Schneider, 2001), pe- 
culiar velocities (Strauss & Willick, 1995) and galaxy rota- 
tion vectors (Lee & Pen, 2000). A map of the density field, 
combined with galaxy surveys, would be a highly valuable 
aid to understanding the formation of galaxies and clus- 
ters of galaxies. Numerical simulations could be performed 
with realizations of initial conditions constrained by our 
knowledge of the density field, allowing object-by-object 
comparison between theory and observation, rather than 
solely statistical comparison. A map could also provide a 
guide to observers who may, for example, wish to search 
for the luminous tracers of the filamentary density struc- 
tures. 

Reconstruction of the density field from galaxy peculiar 
velocities was pioneered by Dekel et al. (1990). The radial 
component of peculiar velocities of galaxies can be deter- 
mined by inferring the distance and then subtracting off 
the Hubble flow contribution to the redshift. Galaxy pecu- 
liar velocity determinations, because they rely on distance 
determinations, are only useful at z < 0.1. 

Peculiar velocities of galaxy clusters determined from 
observations of the Sunyaev-Zeldovich (SZ) effects (Sun- 
yaev & Zeldovich, 1980) do not rely on a distance determi- 
nation. The peculiar velocity signal arises from the clus- 
ter's radial motion with respect to the cosmic microwave 
background (CMB). High-resolution (< 1'), multi-frequency 
observations of galaxy clusters can be used to determine 
the radial velocities of galaxy clusters with an accuracy 
that is independent of the distance to the galaxy cluster^. 

^This independence only fails for clusters at z < 0.2 which suffer 



Errors as low as lOOkms"^ may be achievable (Holder, 
2002). In addition to potential reconstruction, velocity 
measurements will be useful for cosmological parameter 
estimation (Peel & Knox, 2002). 

Here we present a method for rapidly forecasting the 
accuracy of density field reconstruction from peculiar ve- 
locity measurements and apply it to surveys with various 
redshift ranges. We also show how weak lensing obser- 
vations can provide complementary information (Mellier, 
1999; Bartelmann & Schneider, 2001). 

Our analytic method for forecasting results assumes a 
uniform and continuous velocity field map. Even if this 
map were derived from noiseless peculiar velocity mea- 
surements, there would still be an effective noise contribu- 
tion from the fact that most of the mass in the Universe 
is not in galaxy clusters. Fortunately, as we quantify be- 
low, the contribution to the peculiar velocity variance from 
each logarithmic interval in wavenumber k drops as 
for relevant scales, so this noise from under-sampling is 
not overwhelmingly large. 

We consider three different survey types labeled "SDSS" , 
"DEEP/VIRMOS" and "SZ" with redshift ranges 0.2 < 
z < 0.4, 0.7 < z < 1.4 and 0.2 < z < 2 respectively. 
Galaxy cluster redshifts are required to convert SZ mea- 
surements into gravitational potential maps. Thus, the 
redshift ranges of two survey types are subsets of those of 
optical surveys SDSS (Sloan Digital Sky Survey), DEEP 
II (Deep Extragalactic Evolutionary Probe; Davis et al. 
(2001)) and VIRMOS (Visible-Infrared Multi-Object Spec- 
trographs; Le Fevre et al. (2001)). The peculiar velocity 
surveys we imagine could be realized by targeted multi- 
frequency SZ follow-up on galaxy clusters identified in 
those optical surveys. The ambitious, deep "SZ" survey 

significant confusion with primary CMB anisotropy due to their large 
angular sizes. 
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is modelled after proposed SZ surveys which can locate 
galaxy clusters over a large z range. Again, the velocities 
can be determined from multi-frequency follow up (if nec- 
essary). The redshift determinations for this last survey 
type arc left as a currently unsolved challenge. 

Radial peculiar velocities are directly sensitive to ra- 
dial gradients in the gravitational potential. Gravitational 
lensing is caused by tangential gradients in the gravita- 
tional potential; thus, in principle, lensing provides com- 
plementary information for potential reconstruction. How- 
ever, we show below that lensing is actually unlikely to 
add much to the density field reconstruction. The rea- 
son is simple: lensing observations are two-dimensional, 
whereas velocity observations arc three-dimensional. The 
conclusion may be different if one considers lensing of 
many source populations with differing redshift distribu- 
tions (Hu & Keeton, 2002). 

2. FORECASTING ERRORS 

To begin, we assume that we have a continuous and 
uniform map of the velocity field in a box of volume V . 
The map is not noise-free, but has white noise with finite 
weight-per-comoving volume w\ i.e, the variance of the 
error on the average velocity in some sub-volume V\ is 
Oy^ = \l(wV\). For specificity we set our fiducial model 
to be /i = 1 — 0„j = f^A = 0.7 and cg = 1. 

For the moment we ignore the fact we can only make 
measurements along our past light cone and assume the 
map is of the velocity field today. We remove this oversim- 
plification below, but for now this idealization maintains 
homogeneity and allows a completely analytic analysis of 
the potential reconstruction errors in Fourier space. With 
the assumption of a potential flow (Peebles, 1993) 



where 



and 



«;^*A:(io)x(i) 



(1) 
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The noise and signal contributions to the variance of $g 
are diagonal with diagonal entries given by 

N{k) =(^) /w and (4) 
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respectively where A^(fc) = k^Ps{k)/{2w'^) and Ps{k) is 
the matter power spectrum. 

The square of the signal to-noise ratio for the gravita- 
tional potential map is thus 



X 



0.21 
= (100 kms"^' 



^0.1 Mpc" 

uj. =y,^^^.,.n /(64Mpc)^ (6) 

where xq = x{to) = 1.1 and A2(0.1Mpc"^) = 0.8 for our 
fiducial model. We therefore expect measurements with 



S/N > 1 for every mode with 27r/L > fc > 0.1 Mpc"^ 
where = V. Further, there are many of these modes 
for a box of size L; the number in Fourier-space volume 

(Afc)3 is 4 X 10^ ((fc/0.1 Mpc"^)(X/1000 Mpc)(AA;/fc))^ 

We now address our value of Wf above. A fundamen- 
tal limit on w is set by the small fraction of mass in 
galaxy clusters. The comoving number density of galaxy 
clusters with lower mass threshold IO^/i-^Mq at z = 1 
is about one per (64Mpc)"^ volume (Holder et al., 2001) 
which roughly corresponds to the volume of a 40 Mpc ra- 
dius sphere. However, a sphere containing this limiting 
mass in a uniform background only has a comoving ra- 
dius of about R = 9 Mpc. Estimating the velocity field 
averaged over R = 40 Mpc from just one sample of the 
velocity field measured over 9 Mpc introduces an under- 
sampling error with variance (at z = 1) of ((w4o — wg)^) = 
(117kms^^)^. Thus we cannot do better than = 
(117kms-i)2(64Mpc)3 at z = 1. The dashed line in Fig. 1 
shows the redshift dependence of this under-sampling noise 
which is due to evolution in n. According to Holder (2002) 
and Nagai et al. (2002) there is also a fundamental limit 
on how well one can infer the peculiar velocity of a galaxy 
cluster from the (highly non-uniform) velocities of the gas, 
of about lOOkms"^. The result of adding these errors in 
quadrature is the solid line of Fig. 1. Wiener filtering can 
reduce the error, though we do not employ it here. 

The angular resolution, sensitivity and frequency cov- 
erage required to achieve a velocity measurement with a 
given precision have not yet been worked out in detail. 
From prior work (Haehnelt & Tegmark, 1996; Holzapfel 
et al., 1997; Kashlinsky & Atrio-Barandela, 2000; Aghanim 
et al., 2001), which ignores the velocity substructure of the 
gas, we find measurements with errors of 8^K translate 
into (T„ ~ 100 km s"-'^, well-matched to the Holder/Nagai 
limit. A multi-mode detector on a large telescope with 
sensitivity of lOO^Ks^/^ could achieve SjjK on 20 clusters 
in 1 hour, or 16,000 in a month. 

We now consider observations of the radial component 
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Fig. 1. — Inverse weight-per-solid angle, w^^, from under- 
sampling noise (dashed line) and from under-sampling noise plus a 
100 km s^^ error on each galaxy cluster added in quadrature (solid 
line). The dotted line is R(z)/30Mpc where 4/3TrR^ {z)n{z) = 1 
defines R{z) and n{z) is the comoving average number density of 
clusters. 
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of peculiar velocities on our past light cone. The isotropic 
but inhomogeneous geometry suggests the mode decom- 
position 

/•OO 

^(7,r-)=V/ dkk^^im{k)ji{kr)YUj) (7) 
im Jo 

as used by e.g. Stebbins (1996). The line-of-sight radial 
velocity generated by this potential is 

vr{j, r) = ix{r)/Ho V / dkk^^im{k)j'i{kr)Yim{7) 

(8) 

where we have changed our independent variable from t to 
conformal distance, r and j[{y) = d/dy ji{y). Note that 
all the time-dependence on the right-hand side is carried 
by x{r), thus we are discussing reconstruction of the grav- 
itational potential on our past light cone extrapolated by 
linear theory to what it would be today. 
Our discretized model of the data is 

Vi ='^Ay{i;b,l,m)^bim + ni (9) 

blrn 

where n is the error (whose statistical properties were just 
discussed above), Ay is implicitly defined by Eq. 8 and b 
enumerates the k bins: 

him = -TTT / dk^mik). (10) 

i^Kb Jkb-Akb/2 

The minimum-variance estimate of ^um is 
kim= W-i{blm,b'l'm') 

h'l'm'ii' 

X Ay(i-b',l',m')N7.}vi, (11) 
where Nij = {niUj), 

W;^{blm, b'l'm') = {5^bim5^vvm') (12) 
is the noise covariance matrix of ^um given by 
Wy^{blm,b'l'm') = Y^Ay{i^b,l,m)N7^}Al{i'■b',l',m') 

ii' 

(13) 

and 5^blm = ^blm - ^Im- 

With the uniform sampling assumption, W^^ for recon- 
structing ^bim is given by 

W^^iblm, b'l'm') = ^^SwSmm'klkl, 

X / drr'^x^{r)w(r)j'i {kbr)j'i{kb'r). (14) 
The matrix is block-diagonal with zero entries for any el- 
ements with I ^ V ov m ^ m' . As one expects from 
statistical isotropy W^^ does not depend on m. 

The signal matrix is diagonal, though the effect of bin- 
ning the uncorrelated k modes is to reduce the variance by 
Afc: {^bim^Um) = Pq,{k)/^k. In the results section we 
compare signal and noise by plotting k^ /^k/W^^ where 
the k^ factor makes a dimensionless quantity and the A/c 
makes the corresponding signal quantity, k^P^{k) = (47r)^ 
xk^P^{k), independent of Afc. 

The weight matrix for potential reconstruction from a 
weak lensing convergence map with uniform weight per 
solid angle w: 

W^^{blm,b'l'm') = 5w5ram'we{l + If 
xAkbklAkb'k^,Ii{kbrs,e^)Ii{kb>r^,es) (15) 



where 

li{krs,0s) = r — - r)ji{kr)^^p^^W{kr9s) 
Jo rrs a{r]o - r) 

(16) 

(Stebbins, 1996), and WikrOs) = 2J-i,{kre s) / {krO s) (Jain 
& Seljak, 1997) incorporates the fact that we smooth the 
signal over radius 6s = 2.5'. Again, this is block-diagonal 
in I and m. However, the structure in the radial wave- 
number k is highly degenerate. For each l,m sub-matrix, 
the fc, k' dependence can be written as an outer product 
of a single vector. Such matrices have only one non zero 
eigenvalue, a manifestation of the inability of lensing mea- 
surements to distinguish between different radial modes. 

3. RESULTS 

In Figure 2 we compare signal and noise for AA; = k for 

several values of / assuming w{z) as shown in Fig. 1. One 
can sec that S/N > 1 measurements arc possible on large 
spatial and angular scales. The sharp low k cutoffs are an 
effect of the finite size of the surveys: fcmin = l/L where L 
is the radial depth of the survey since ji{kr) ~ for A; < 
l/r. The velocity results can be understood quantitatively 
by multiplying S/N from Eq. 6 by the number of modes 
in a bin of width AA:: AAL/(27r). 

Now we turn to the weak lensing panel of Fig. 2. We 
assume a noise level of w; = 9 x 10^ which can be achieved 
with a galaxy number density of n = 30 gal. arcmin"^ 
and intinsic galaxy ellipticity of (Tc = 0.2 (Kaiser, 1998; 
van Waerbcke et al., 1999). We further assume a negli- 
gible width to the distribution of these galaxies centered 
at Zs = 1. The results can be understood analytically 
from Eq. 15 and the approximation (Stebbins, 1996) /; = 
^/'K/{2l){l/l - l/{krs)) for krg > I and /; = otherwise. 
Unlike with velocity the higher / values are reconstructed 
with less noise due to the sensitivity of the convergence to 
tangential gravitational potential variations. 

Above we have plotted kf /W^^ {blm, blm) as an indica- 
tor of the variance expected on a measurement of kf^bim- 
However, the variance of the error on this mode is actu- 
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Fig. 2. — Dimensionless diagonal noise variance {Akk^ /W-^^) for 
binned modes of the gravitational potential with bin widths AA; = k 
for several values of I, for our three velocity surveys (X = v) and 
the lensing survey {X = k). The expected signal variance, fc^P,j(fc) 
is also shown for comparison. 
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ally given by kfW~\{blm,blm). This inverse is not well- 
defined because of eigenmodes with zero eigevnalues. In 
Fig. 3 wc plot the eigenvalues of the signal-to-noise matrix 

W'£p^/AkWl^' (e.g., Bond (1995)). 

The total weight from a combined weak lensing and pe- 
culiar velocity survey is simply given by adding the indi- 
vidual weights. The eigenvalue spectrum for this combina- 
tion is identical to the peculiar velocity survey eigenvalue; 
spectrum, except for the increased value of one eigenvalue 
per l,m mode. 

4. DISCUSSION 

The largest / value of modes with S/N > 1, /max, entails 
a minimum survey size of 6*111111 = 5°(80/Zniax)- From Fig. 3 
wc see ^max — 80 and therefore 9 

min 

5° for the two 

deeper surveys and Imax — 35 and therefore ^min — 10° for 
SDSS. 

There are many advantages to the shallower SDSS sur- 
vey. The spatial number density of clusters is higher, re- 
ducing the under-sampling noise. Despite the increased 
value of ^mim the clusters in solid angle 0^^^ are actually 
fewer for the shallow survey (358) than for the deeper sur- 
veys (500 and 1000). This is an advantage since at fixed 
detector sensitivity the requisite observing time is propor- 
tional to the number of clusters, not the survey solid angle. 
Nearer clusters are easier to observe in the optical and in 
X-rays, so redshift and temperature determinations will 
take less time. Finally, a nearby map would be easier to 
compare to other tracers of the density field than a more 
distant map would be. 

Increasing the survey solid angle 9"^ not only reduces 
^min, but also iucrcascs spectral resolution (i.e., decreases 
SI = 2tt/9 = 12(30°/^)) and increases the number of an- 
gular modes with S/N > 1: 



l=l„ 



9 



30' 



35 



TT. (17) 



Thus for a 61 = 30° survey with 0.2 < z < 0.4 the total 
number of modes with S/N > 1 is roughly iVg x 15 ~ 
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Fig. 3. — Signal-to-noisc eigenvalue spectra for individual angular 
modes. The boxes show the single non-zero eigenvalues of the weak 
lensing survey. The SDSS panel includes additional curves for I = 10 
(top curve) and I = 50. 



360. Covering the entire tt steradians of the SDSS survey 
would require measuring 36,000 clusters and would deliver 
~ 4, 500 modes. 

5. CONCLUSIONS 

We have shown that galaxy cluster peculiar velocities 
can be used to make high signal-to-noise maps of the 
gravitational potential, and therefore matter density, on 
very large scales. A limiting source of uncertainty in these 
maps is what wc have called the under sampling noise due 
to the fact that only a small fraction of mass is in clusters. 
The under-sampling noise can be reduced by lowering the 
mass threshold but the lower the mass of the cluster, the 
harder the measurement of its peculiar velocity. Much 
more work is needed to understand the demand on experi- 
mental resources and to optimize observing strategies. We 
have made several arguments that favor shallower surveys 
over deeper ones. 

The $ and density maps will have numerous applica- 
tions. Cross-correlations with various tracers of the den- 
sity field (red galaxies, blue galaxies, infrared galaxies, 
lensing of galaxies, lensing of the CMB, integrated Sachs- 
Wolfe effect on CMB, ...) will be of great interest. The 
greatest drawback of the maps is their inability to probe 
scales below the typical separation between clusters. 

We thank G. Holder for his encouragement and for shar- 
ing the n{z) calculation used in Holder et al. (2001) and L. 
van Waerbeke and S. Colombi for several useful conversa- 
tions. O.D. acknowledges financial support from "Societe 
de Secours des Amis des Sciences" . 
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